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Dust Scattering Models for T Tauri Disks and Envelopes 

With the Monte Carlo code developed by B. Whitney and L. Hartmann, we have 
computed a series of models of scattering in disks around young stellar objects. The 
code calculates scattering by dust, including polarization, in arbitrary geometries. By 
computing model images, we found that disks, by themselves, around young stellar 
objects would be very difficult to detect with present-day imaging techniques. The 
canonical geometrically thin accretion disk illuminated by a central T Tauri star scatters 
very little light at distances of many AU from the star. An optically thick disk which has a 
“flaring” surface may be observable, especially if viewed nearly edge-on so that the stellar 
source becomes occulted. Such a surface provides more solid angle to intercept (and 
scatter) stellar light. Example images are shown in Figure 1 of a flared disk illuminated 
by a star at the center. The disk is 500 AU in radius and the scale height of the disk 
increases with radius as h ~ r 9 ^ 8 . The contours show the light distribution at V (5500 
A), normalized to the peak flux in a pixel of the image. Each contour level increases by 
a factor of 2. The minimum contours for each image are, going in increasing /x, 0.01, 

1 x 10~ 4 , 5 x 10~ 5 , and 5 x 10" 5 . The straight lines indicate polarization vectors, with the 
position angle shown by the orientation of the lines and the length of the lines showing 
the magnitude of the polarization. The total polarization is shown with each image. 
The position angle is oriented perpendicular to the disk plane. Only the image in which 
the central source is occulted (/x = 0.1) shows a high amount of relative flux, and high 
polarization. The other images have such high contrast between the stellar flux and 
diffuse flux that measurement would be difficult. 

In comparing these images to observations of young stellar objects which show diffuse 
structure, we find little resemblance. Objects such as HL Tau and L1551 IRS 5 show 
a bowl-like structure at infrared wavelengths, which has been attributed to a flared 
disk by several authors (Monin et al . 1989, Strom et al. 1985). These objects also 
have high polarization with position angle oriented parallel to the presumed disk plane 
(perpendicular to the jet directions). The high polarization indicates that the central 
source is obscured, and the position angle indicates that significant material lies above 
the disk plane. A flared disk system will only give high polarization when viewed edge-on, 
and the position angle is always oriented perpendicular to the disk plane. 

This suggests that an envelope, perhaps the remnant infalling envelope, must be 
present to scatter more stellar light than a disk can, and obscure the star at many 
inclinations. We have computed a grid of models of scattering in a disk+envelope 
system. We use the envelope density structure of a rotationally flattened cloud in free fall 
(Terebey, Shu and Cassen 1984, Cassen and Moosman 1981). The addition of angular 
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momentum to the infall solution causes material to fall onto a disk thereby decreasing 
the optical depth through the envelope to the star. This allows scattered light to be 
seen even in embedded protostars assuming the usual values of infall rate and rotation 
rates for the Taurus molecular cloud. Wind-blown cavities allow stellar light to escape 
more easily and scatter off the walls of the envelope. This is illustrated in Figure 2, 
which shows an object surrounding by an envelope with bipolar holes. The hole shape 
was chosen to be a streamline of the material which falls in from an angle with cosine 
Ho — 0.87. We assume a wind has blown out a hole within this angle and particles 
continue to fall in on streamlines beyond this angle. At large distances the streamline 
shape is radial, as can be seen in Figure 2. The images are shown at an inclination of 
H = 0.6. The left panels are models run with forward-scattering high-albedo dust, typical 
properties at visual wavelengths. The right panels have low-albedo, Rayleigh scattering 
dust more likely in the infrared. The opacities vary by a factor of 10 going from top to 
bottom, as shown by the optical depth variation. The optical depth through the envelope 
is proportional to K/(250cm* /gm) M/(2 x lO^N^yr" 1 ) (^/(lOOAU))- 1 / 2 , where r c is 
the centrifugal radius determined by the angular momentum of the original cloud at the 
furthest radius from which the material collapsed. Thus the sequence in images from top 
to bottom can be thought of as a sequence in wavelength (opacity) going from optical to 
IR, infall rate, disk size, or combinations of these. 

The top panels have an optical depth of 9.1 through the envelope yet a significant 
amount of scattered light can be seen. The polarization maps resemble observations 
of many young stellar objects (e.g. Gledhill and Scarrott 1989; Draper, Warren-Smith 
and Scarrott 1985; Scarrott, Draper and Warren-Smith 1989), especially the elliptical 
pattern of the vectors at the higher optical depths. The total system polarization is 
aligned parallel to the disk at high opacity and perpendicular at low. This position 
angle rotation is a feature observed in some objects such as T Tau. Note that the 
polarization can be large even in very optically thick envelopes, if it has wind-blown 
holes. A typical dust albedo of u> < 0.5 insures that most photons are singly scattered, 
and a hole allows light to scatter mostly in the polar regions. Such an asymmetric 
geometry for scattering leads to very high polarization. Thus, large polarization should 
not be interpreted as due to single scattering in optically thin envelopes. In fact, at low 
optical depths, the envelope becomes more symmetric, the star less obscured, and the 
total system polarization actually decreases (lower panel in Figure 2). 

The longer wavelength images appear bowlshaped. This is not due to scattering in 
the disk but scattering in the envelope. The objects that typically show these “disk”-like 
structures are usually the same ones with far-IR excesses, fit by accretion disk models with 
temperature-radius relations which are too flat to be explained with normal accretion. 
The spectral energy distributions of these objects should perhaps be re-interpreted in the 
context of envelopes instead of accretion disks. 
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FU Orionis Disk Winds 


The discovery of powerful bipolar outflows from pre-main sequence objects indicates 
that mass loss plays an important role in star formation (Lada 1985; Shu, Lizano, &; 
Adams 1987; Edwards, Ray, & Mundt 1991; and references therein). Bipolar outflows 
appear to be quite common in the earliest stages of star formation (Terebey, Vogel, k 
Myers 1989) and can easily disperse the material surrounding the central source. In 
certain cases, long and remarkably thin optical “jets” of material emanate from young 
objects (Reipurth 1989). The kinetic energy carried by these jets often amounts to ^ 
10% of the total system luminosity (Lada 1985). 

The most popular theories for producing bipolar flows or jets from young stars tap 
the energy liberated as material accretes through a circumstellar disk onto the central 
star (Adams, Lada, k Shu 1987; Kenyon k Hartmann 1987; Bertout, Basri, k Bouvier 
1988; Basri and Bertout 1989). In these theories, high velocity gas is ejected either from 
a rapidly-rotating central star (Hartmann k MacGregor 1982; Shu et al. 1988; Ruden 
et al. 1990) or from an accretion disk (Pudritz k Norman 1983; Konigl 1989; Pringle 
1989; Pelletier k Pudritz 1991). In both types of models, material is flung outward along 
strong magnetic field lines anchored in the star or in the disk. The rotating magnetic 
field ultimately collimates the flow into bipolar jets (e.g., Blandford and Payne 1982). 

The source of the mass loss from young stars is controversial, because many bipolar 
flow sources are heavily reddened and therefore not directly observable. However, optical 
observations can be employed to study the wind acceleration region in several young FU 
Orionis variables that have powerful outflows and are not strongly reddened (Reipurth 
1989, 1990; Hartmann et al 1991). FU Orionis objects are thought to be protostellar 
disks undergoing rapid mass accretion at rates of M^cc ^ 10 -4 Mq yr -1 (Hartmann k 
Kenyon 1985; Hartmann et al. 1991), and thus provide good laboratories for studying the 
connection between protostellar disk accretion and mass loss. The prototype of the class, 
FU Ori, has a strong wind (Herbig 1966; Bastian k Mundt 1985), with a mass loss rate 
~ lO _ 5 M 0 yr _1 (Croswell et al. 1987), and a terminal wind velocity v ~ — 300 km s -1 
detected in the strong Ha line. 

If FU Ori has a disk wind, spectral lines of different strengths should display different 
profiles. Strong lines, formed in the outermost disk atmosphere, should exhibit strongly 
blueshifted absorption components. Conversely, weak fines formed in deeper, denser disk 
layers should exhibit profiles in which rotation dominates, since by mass conservation 
the expansion velocities will be small in dense regions. The fine profiles should vary 
continuously as the fine opacity increases, so that fines of intermediate strength should 
exhibit intermediate (small) blueshifts. Observations provide qualitative support for 
this picture; strong Fe II lines exhibit “shell” absorption at velocities of ~ —80 to 
— 50kms -1 (Herbig 1966; Petrov k Herbig 1992), while weak fines are symmetric and 
double-peaked. In addition, moderately strong photospheric fines in the blue spectral 
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region are asymmetric, and show blueshifts of — lOkms 1 (Hartmann & Kenyon 1985, 
1987a, b; Kenyon et al. 1988 = KHH; Welty et al. 1991). Although this behavior is 
consistent with the disk wind hypothesis, no quantitative model for the lines has yet 
been constructed. 


We have submitted a manuscript entitled “Mass Loss from Protostellar Accretion 
Disks. I. The Accelerating Wind of FU Orionis” by Calvet, Hartmann, and Kenyon. 
In this paper we present evidence that the wind of the pre-main sequence object FU 
Orionis arises from the surface of the luminous protostellar accretion disk. A disk wind 
model calculated assuming radiative equilibrium explains the differential behavior of the 
observed asymmetric absorption line profiles. The model predicts that strong lines should 
be asymmetric and blueshifted, while weak lines should be symmetric and double-peaked 
due to disk rotation, in agreement with observations. 


We considered a viscous, geometrically thin, optically thick accretion disk with 
constant mass accretion rate surrounding a star with mass M and radius /?,. The 
reference height H measured in the direction outward from the disk midplane (z) varies 
with cylindrical radius r a s H = H 0 (r/Ri) 9/8 . We take H 0 = 0.1/2,, which is consistent 
with some theoretical models of FU Ori objects (Clarke et al. 1990). 

Our radiative transfer solutions assume that spectral lines form close to the (geo- 
metrically thin) disk, so we neglect radial energy transport and adopt a plane- parallel 
geometry. We first calculate a plane-parallel model atmosphere at each disk annulus 
assuming hydrostatic equilibrium, as discussed in detail by Calvet et al. (1991b). This 
model neglects viscous dissipation within the photosphere, so radiative equilibrium holds 
and viscous energy generation in deeper disk layers determines the surface radiative 
flux. The constant radiative flux through the optically thick disk atmosphere can be 
characterized by an effective temperature 


leff = 


' 3GM M» 

47T<7T 3 


nr)] 
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( 1 ) 


(Shakura & Sunyaev 1973; Lynden-Bell & Pringle 1974). The vertical temperature 
structure is calculated using the grey atmosphere approximation in the Eddington 
approximation for a Rosseland mean optical depth t. To calculate the structure of 
the atmosphere at each radius r we assume that the reference height H corresponds to 
the level r = 2/3, and solve by iteration the system of equations 

dp _ gjz) 

dr XRoss ' 
dz 1 

df XRosi P 

where the height-dependent gravity is given by 

GM ( H + z) 
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(3) 


9{z) = 


[r 2 -f (z + /f) 2 ] 3 / 2 


(4) 


4 



Here XRo*» is the Rosseland mean opacity, p is the density, and z is the height of the 
atmosphere above the r = 2/3 level. Atomic and molecular populations are calculated 

assuming LTE. 

Having calculated a radiative equilibrium atmosphere assuming hydrostatic equi- 
librium throughout, we then modify the structure by introducing the disk wind. For 
simplicity we assume that material flows away from the disk midplane m the vertical 
(z) direction at a velocity v. We further assume that the mass loss per unit area pv at 
cylindrical disk radius r varies as 

where M 0 and cq are constant input parameters. At each r we define a level z, from the 

relation . 

(pv) r (z s ) - P, c„ v 6 ) 

where p, is the density given by the hydrostatic equilibrium solution and c, is the 
adiabatic sound speed (7 = 5/3). We assume that the flow is smoothly accelerating 
outwards; therefore the flow is subsonic for z < z, and is supersonic for z > z,. In the 
subsonic region, below z„ we assume that the density structure is given by hydrostatic 
equilibrium, which should be true asymptotically proceeding toward the disk midplane. 
In the supersonic region, above z 4 , we adopt a velocity law u(r, z) and compute the 
density from mass conservation for steady flow. 

The velocity structure of a rotating hydromagnetic wind depends in detail on several 
parameters including the magnetic field strength and structure that are not currently 
predictable from first principles (Blandford k Payne 1982; Lovelace et al 1987). The 
flow velocities of the spherically self-similar models of Blandford and Payne (1982) scale 
as V(z/r) oc V Kep {r 0 ), where V Kej> is the Keplerian velocity where the field line enters 
the disk at r = r 0 . We therefore adopt a simple parameterization of the vertical velocity 
structure and assume that the wind velocity increases linearly with distance above z„: 

u(r, z) = T V Kep (r ) x (z - z.)/r + c s {z ,) , * > • ( 7 ) 

The density in the supersonic region is then given by 

f x (Hr (8) 

^ (r ’ z) "v (1 + 1.58 (z + H)fr)' 

with « determined by equation (7), and z m0I is an arbitrary cutoff chosen so that 
v(z max ) = 100 kms' 1 . (The factor (1 + 1.58(z + H)/r) represents a departure from mass 
conservation in cylindrical geometry, and is motivated by the divergence of field lines 
in the standard model of Blandford k Payne (1982). The difference from cylindrical 
mass conservation only becomes important when z w r, where our approximation of 
plane-parallel geometry breaks down anyway.) 
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With fixed p(z), the new Rosseland optical depth and temperature scales are deter- 
mined at each 2 by integrating inward from z max , 

t(z) = T{z - A 2) + x(p(z), T(z))Az . ( 9 ) 

Using the grey Eddington expression for the temperature, 

W = 3/4 T!„ (r(2) + 2/3), (10) 

equation (9) can be solved as an implicit equation for t{z). Finally, the pressure is 
calculated from the equation of state, interpolating at given density and temperature in 
p re-calculated tables of mean molecular weight. 

The introduction of the wind, with finite continuum optical depth, increases the 
total optical depth in the subsonic region and thus increases the temperature for 2 < 2,. 
Consequently, the hydrostatic equilibrium solution for this region is no longer self- 
consistent; but for most cases of interest, the change in the temperature is unimportant. 
Similarly, the change in the optical depth scale due to addition of the wind also changes 
the effective gravity of the photosphere, since <7(2) is height-dependent, but our results 
are only sensitive to changes in gravity of an order of magnitude. 

The radiative transfer models depend upon the adopted disk and wind parameters. 
Observational estimates of the central mass, inner disk radius, and accretion rate in FU 
Ori are dependent upon the inclination angle, which is only poorly constrained (KHH). 
We adopt the accretion disk parameters used in our treatment of the near-infrared 
spectrum of FU Ori (Calvet et al. 1991a): M = 0.49 M 0 , R, = 4.42 Ro, = 
1.59 x 10 M e yr *, except that we take i = 35° rather than i = 40° to provide a slightly 
better match to rotational velocities. These parameters fix the effective temperature 

(equation [1]), surface gravity, and the projected rotational velocity of the disk as a 
function of radius. 

The final spectra were computed using a spectral line list kindly provided by R. 
Kurucz. The gf values were adjusted to provide agreement between the spectrum of a 
G supergiant and a non-rotating model atmosphere of appropriate effective temperature 
and gravity. (We used a G supergiant to provide the best match to the average optical 
spectrum of FU Ori; see KHH). 

In the context of the present simple model, the wind properties are fixed by the 
parameters M 0 , a u and I\ The mass loss parameter M 0 is constrained by the results 
of Croswell et al. (1987), who estimated a total mass loss rate for the wind M w ~ 
10 M 0 yr for FU Ori from modelling the Ha and Na I resonance lines. In our final 
models, we adopt M 0 ~ 1 x 10 -5 M 0 yr -1 . The other two wind parameters are not as well 
constrained by other considerations. We adopt a x = 2, as suggested by the self-similar 
wind solutions of Blandford & Payne (1982) For this choice, the total mass loss rate from 
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the disk at radii interior to r is 


M w (r) = M„ln(r/Ri). ( n ) 

Finally, T was adjusted to produce the best agreement with observations. We found that 
T = 1 results in acceptable line profiles, and is roughly consistent with the velocity law 
in the standard model of Blandford & Payne (1982). 

The wind temperatures of the radiative equilibrium models are in reasonable agree- 
ment with the results of Croswell et al. , who estimated T w ~ 6000 K for FU Ori; for 
comparison, at r = 2R,, T w asymptotically approaches 5420 K. At large distances above 
the disk, where gas densities are low, there may be departures from radiative equilibrium 
caused by adiabatic expansion and possible chromospheric heating. We ignore these 
effects for our initial exploration. 

At the innermost disk radius we take to have a wind, r = 136 Ri, a mass loss rate 
of M 0 = 10 -5 M© yr -1 has a profound effect on the atmospheric structure. Because the 
continuum opacity depends strongly on temperature, even the relatively modest increase 
in T efJ from 6440 K at 2 R, to 7190 K at 1.36 R, causes the photosphere, z(t ~ 2/3), to 
move to considerably higher vertical distances. More importantly, the expansion velocity 
at the photosphere exceeds lOkms" 1 . Thus, at this mass loss rate, and with T = 1, all 
photospheric lines formed at this radius will be appreciably blueshifted. 

Roughly 60% of the optical spectrum arises from disk annuli between 1.5 R, and 3 R, 
(KHH), which can be represented adequately by our results for r = 2R;. The vertical 
structure at this annular distance begins to depart considerably from the purely static 
atmosphere when t Ro33 <> 2xl0~ 2 . The wind velocity, £ 30 kms" 1 , is already a significant 
fraction of the disk’s projected rotational velocity of ~ 6010 ns- 1 , so lines formed at or 
above this level will exhibit profiles strongly affected by expansion. 

To illustrate the effects of the wind, we consider three spectral .lines which have 
qualitatively different behavior. The Li I 6707 A line is always double-peaked in spectra 
of FU Ori (Hartmann et al. 1991; Petrov & Herbig 1992; Welty et al. 1992). The strong 
Fe I 4957 A line is generally not doubled, but roughly parabolic in shape, with little 
or no velocity shift (cf. Petrov & Herbig 1992; Welty et al. 1992). The Fe II 5018 A 
line shows appreciable blueshifted absorption in a shell feature. We wish to explore 
whether a wind with a reasonable mass loss rate can explain the differences in behavior 

among these lines. 

The emergent spectrum is the sum of spectra calculated for different annuli with 
differing rotational velocities and wind velocities. Figure 3 compares the resultant line 
profiles with observations of FU Ori obtained with the KPNO 4m echelle spectrograph 
and TI CCD in December of 1988. Further details of the observations and data reduction 
can be found in Hartmann & Stauffer (1989) and Hartigan et al. (1989). 
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The model predicts that the Li I line in FU Ori should be double-peaked, as observed, 
because it forms in a deep layer where the dominant motion is Keplerian rotation. We 
increased the Li abundance to log Nn = 3.6 (on a scale where log Nfj = 12), a factor 
of two or more increase over typical pre-main sequence abundance estimates, to provide 
a better match with observations. The calculated equivalent width is still too small, 
but the Li I line is very temperature- sensitive, and a modest change in the disk effective 
temperature might substantially improve the agreement with observation. 

The Fe II A 5018 line lies in the blue spectral region, where line blending is extremely 
important. The effect of the wind is clearly demonstrated by comparison with the model 
with no mass loss. The wind model exhibits an absorption “dip”, as observed. In our 
model, this absorption dip is not the result of a “shell” of material at nearly constant 
velocity (e.g., Herbig 1966), but instead arises from the convolution of rotation with the 
wind line profile, as described in the previous section. 

We note that agreement between theory and observation requires strong Fe II 
absorption extending to blueshifts of u mox = 100 kms" 1 . At this velocity, with T = 1, 
z ~ r , the assumption of plane-parallel geometry breaks down. However, our model 
clearly predicts that the Fe II line must be very asymmetric, as observed. 

The results for the Fe I 4957 A “line” are complicated to interpret because this 
feature is actually a blend of several lines. Two strong Fe I lines from the same multiplet 
are present; the 4957.60 A line dominates the blend, but the 4957.30 A line cannot be 
neglected. The blend causes the profile of the 4957 A feature to appear only slightly 
doubled, even in the absence of any wind (Figure 3). 

Our final model for the 4957 A feature, shown in Figure 3 for a mass loss rate of 
1 x 10 M@yr 1 , produces a line blend that is not doubled, as observed, but has a larger 
equivalent width than indicated by the data. We feel that the agreement between theory 
and observation is reasonable, given the difficulties of modelling the blended spectrum. 
In particular, if the model continuum were depressed by approximately 10 % by weak line 
blanketing, the agreement between theory and observations would become much better. 

We emphasize that our model results for the 4957 A feature are the product of 
several adjustable parameters, such as log gf values, turbulent velocities, and blends with 
unidentified features, in addition to the basic model paramters M 0 and T. Thus, our 
calculations for this region are not definitive. However, the models demonstrate that for 
plausible wind parameters, the shape of the 4957 A feature can be affected by mass loss 
in such a way as to suppress the double-peaked structure, as observed. 

Petrov &£ Herbig (1992) studied the profiles of about 50 reasonably unblended 
(mostly neutral metallic) absorption lines in the blue spectrum of FU Ori. Many of 
these lines exhibited two absorption components. Petrov & Herbig showed that the 
velocity difference between the two absorption components decreases with decreasing 
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residual line intensity (increasing line strength). In addition, Petrov k Herbig found 
that the long-wavelength absorption component of becomes progressively blues 1 te 
with increasing line strength, while the short-wavelength absorption component stays 
roughly at the same velocity (their Figure 7). 

The wind model qualitatively explains these effects. To show this without laboriously 
fitting every line blend, we constructed a sequence of line profiles for Mo = 10 5 M© yr , 
using the wavelength and excitation parameters of the Fe I 4957.6 A line, but using 
different gf values. This should provide a qualitative treatment of the profiles of neutra 
metallic lines with excitation potentials ~ 2.8 eV (typical of many lines in this spectral 

region). 

In Figure 4 we display this sequence of Fe I profiles for differing gf values. The 
velocity splitting between the two components decreases as the line strength increases. 
In addition, the blue absorption minimum remains at nearly the same velocity, 40 to 
— 35kms _1 , while the red absorption minimum moves considerably closer to line center. 
Both of these effects are qualitatively similar to the results reported by Petrov k Herbig 
(1992). When the line is very strong, only one absorption minimum is observed, also in 
agreement with the results of Petrov k Herbig. It is difficult to make these comparisons 
more quantitative, because small details of the line profiles depend sensitively on line 
blending, but it appears that the models reproduce the qualitative behavior of the 

observations. 

Spectral observations at infrared wavelengths sample outer, cooler regions of the 
disk, compared with optical observations (KHH). Thus, mass-loss signatures in the 
infrared spectrum provide information on the variation of mass outflow with radius. 
Using cross-correlation analysis, Hartmann k Kenyon (1987a, b) and Kenyon k Hartmann 
(1989) suggested that the v" - v' = 2 - 0 first overtone absorption band of CO at 2.3 
am might be affected by mass loss. We calculated the 2.3 fim CO. spectrum, using 
line parameters described in Calvet et al. (1991a). We consider the first overtone 
v " -v' = 2-0 and v" - v' = 3 - 1 bands of CO. 

The introduction of a wind with M 0 - 10~ 5 M© yr" 1 produces a small asymmetry in 
the v" - v' = 2 - 0 band. The 4m FTS spectrum of FU Ori is too noisy for individual 
lines to be compared with the model. We therefore analyzed the average CO line profiles 
using cross-correlation analysis (cf. Hartmann k Kenyon 1987a, b; KHH). We used model 
calculations for no wind, * = 0° as the template spectrum for both FU On and the wind 

model. 

The cross-correlation results are shown in Figure 5. The model predicts a modest but 
significant asymmetry of the v"-v' = 2-0 cross-correlation, indicating extra blueshifted 
absorption. This asymmetry is in reasonable agreement with the cross-correlation of the 
observed spectrum. The agreement with observations of the v" - v - 2 0 ban 


9 



may be slightly better for M 0 = 10-*M @ yr-\ but it is not clear whether the limited 
signal- to-noise values of the observations can support such a distinction. 

The cross-correlation of the v" - v' = 3 - 1 region in FU Ori is almost perfectly 
symmetric; the model cross-correlations are only slightly asymmetric, with a larger 
asymmetry at higher mass loss rates. As Kenyon k Hartmann (1989) suggested, the 
first excited vibrational state of CO has a lower population than the ground state in the 
cool wind, and so the wind is more transparent in the v" — v' = 3 — 1 band. 

We find that mass loss rates M 0 < 10~ 6 M 0 yr _1 do not provide any measurable line 
asymmetry in the v" — v' = 2-0 band. The reason for the discrepancy between our 
results and the estimate M w ~ 10 _7 M®yr _1 of Kenyon k Hartmann (1989) is that the 

latter failed to take the distribution of CO molecules over many rotational levels into 
account. 

The model calculations for the near-infrared CO lines are consistent with the wind 
parameters M 0 = 5 x lO' 6 - 1 x 10' 5 M 0 yr' 1 , a, = 2, T = 1, that is, roughly the same 
parameters as were used to fit the optical lines. However, it should be emphasized that 
the observational indications of mass loss in the CO first overtone lines are marginal, and 
that we cannot rule out the possibility that the asymmetry of the v" - v' = 2 - 0 band 
arises from absorption in a distant, cold, low-density circumstellar shell, rather than a 
disk wind (see discussion in Kenyon k Hartmann 1989). 

Our results suggest that the outflow of FU Ori originates over a substantial area of the 
inner accretion disk, but not over many decades in radius as suggested by the Pudritz k 
Norman and Konigl models. If dust is not destroyed at temperatures cooler than < 10 3 K, 
a massive wind produces unacceptably large extinction and infrared excess. This implies 
that the mass loss rate falls off much faster with radius than in the self-similar model 
beyond about 12 /fc. If the blueshifted absorption in the v"-u' = 2- OCO lines arises 
in a distant circumstellar shell rather than a wind, significant mass loss is limited to a 
region of several R, near the central star. 

Estimates of mass loss rates depend upon the choice of velocity gradient. We have 
not been able to reproduce the FU Ori observations with velocity gradients T ~ unity 
for M 0 = 10 6 M 0 yr -1 . Conversely, the mass loss from the optical region of the disk 
cannot be as large as lO _4 M 0 yr -1 , because then the wind would be optically thick 
in the continuum. In this case ALL of the absorption lines would be blueshifted by 
tens of kms- 1 , whereas the FU Ori objects are all within a few kms" 1 of the neighboring 
interstellar medium (Herbig 1966, 1977; Hartmann k Kenyon 1987a, b). Thus, our results 
reinforce the estimate of M w ~ 10" 5 M 0 yr' 1 found by Croswell et al. (1987) from a study 
of the Ha and Na I line profiles. 

The ratio of the kinetic energy loss in the wind L w to the accretion energy is 

Tu;/L acc = (M W V 3 /2L acc ) . (14) 
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Our estimates of M w ~ 10- 5 M© yr' 1 and a terminal wind velocity V w ~ 300 km s 1 imply 
L w ~ 0.3 L ace . This result is an order of magnitude estimate, but it does suggest that the 
wind may be an important component of the energy budget of the inner disk. 

The ratio of wind angular momentum loss to the angular momentum transfer required 
for steady accretion at M«x is more difficult to estimate. Averaged over the inner disk, 
this ratio is 

JjJacc ~ (Mw/M.cc)^/^! 2 ), (15) 

where is the Alfven radius of the wind and ft* is an average Keplerian angular 
velocity in the inner disk. If we take Qr A ~ V w and ft ~ (Pudritz k Norman 1986), 
M w ~ lO^Maco and ft* r ~ 100 km s' 1 , the results suggest that the wind may carry 
away an appreciable fraction of the inner disk angular momentum. However, the limits 
on mass loss at large radii discussed above indicate that the wind is very unlikely to carry 
away most of the angular momentum from the outer disk. 

If the self-similar scaling holds to the smallest radii, we find that the wind from the 
innermost disk regions is optically thick in the continuum. This material would help to 
absorb or “hide” any “boundary layer” emission from the inner edge of the disk. In simple 
disk models, the boundary layer is the narrow disk region where material in Keplerian 
rotation loses enough energy to come to rest on the slowly- rotating star (Lynden-Bell and 
Pringle 1974). It is interesting to speculate whether the absorption of boundary layer 
emission in a optically- thick wind over an area ~ Rf could account for the non-detection 
of boundary layer radiation from FU Ori objects (Kenyon et al. 1990). 
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Figure Captions 

Figure 1 - Images of a flared disk illuminated by a central source, at different inclinations 
p = cos 0 (p = 0 for edge-on). 

Figure 2 - Scattered light images of a disk+envelope system with polar holes illuminated 
by a central source, at an inclination of p = 0.6. The contours increase by factors of 2. 
The minimum contours for each image are, going from left to right and top to bottom, 
0.1, 2 x 10- 2 , 1 x 10- 3 , 5 x 10- 4 , 5 x 10" 4 , 1 x 10" 4 . 

Figure 3 - The profiles of three typical lines, Li I A 6707, Fe I A 4957, and Fe II A 5018, 
over the entire disk, with M 0 = 1 x 10' 5 M 0 yr" 1 and i = 35°, compared with observations 
of FU Ori. 

Figure 4 - A sequence of disk model fine profiles for M 0 = 1O~ 5 M0 yr -1 , for the wavelength 
and excitation parameters of the Fe I 4957.6 Aline, but a range of gf values to illustrate 
the dependence of profile shape on line strength (see text). 

Figure 5 - Cross-correlations of selected CO first overtone band regions. The calculated 
model spectrum without wind and rotation was used as a template both for the wind 
models and for FU Ori. The top panel is the correlation of the spectral region between 
2.291 pm and 2.310 pm, which includes much of the v" — v' = 2 — 0 band away from the 
bandhead; the lower panel is the correlation of the region between 2.321 pm and 2.349 
pm for the v” — v' = 3 — 1 band. 
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